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1.  INTRODUCTION  AND  SUMMARY. 

Let  (N^.Ng)  be  a  stationary  bivariate  point  process  on  R.  This  article  is 

concerned  with  statistical  methods  for  discovering  and  quantifying  an  association 

between  the  two  processes  from  a  realization  A,  <  A- <  . . .  <  A  ,  B,  <  B.  <...  <B  over 

12  n  '  1  2  n„ 

A  B 

a  long  period  of  time  T.  The  paper  is  motivated  by  certain  problems  that  arise  in 
neurophysiology,  which  are  very  briefly  described  as  follows  (for  further  details  see 
e.g.  Bryant,  Ruiz  Marcos,  and  Segundo,  1973). 

Two  neurons.  A  and  B,  are  monitored  over  a  period  of  time  T  during  which  each 
neuron  fires  a  sequence  of  impulses.  The  problem  is  to  determine  whether  or  not  the 
impulse  times  are  associated.  An  association  between  NA  and  Ng  may  be  construed  as 
evidence  that  either  the  two  neurons  are  communicating,  or  that  they  both  share  input 
from  a  third  source. 

Another  problem  arises  in  certain  neurophysiological  studies  of  learning  and 

memory.  An  animal  is  to  be  taught  (trained)  to  perform  a  certain  task.  Now  consider 

two  "connected”  neurons,  A  and  B,  which  are  essential  in  the  performance  of  this  task. 

Record  the  impulse  times  during  a  period  before  the  learning  experience,  obtaining  a 
Bef  Bef 

realization  of  (NA  *,Ng  ),  and  during  a  period  of  time  well  after  the  learning 

experience,  obtaining  a  realization  of  (NA^t* jNg^1’) .  The  processes  N®6^'  and 

Ng6^'  may  be  dependent.  The  problem  is  to  determine  whether  or  not  this  dependence 

Aft  Aft 

is  "stronger"  for  the  processes  NA  ’  and  Ng^  '  .  A  neurophysiologist  may  consider 
a  change  in  the  strength  of  the  dependence  as  evidence  that  learning  has  taken  place. 

The  two  problems  have  very  different  statistical  character.  Let  S  be  a  statistic 
that  "measures"  the  dependence  between  two  point  processes.  The  first  problem  is 
one  of  testing  the  hypothesis  that  NA  and  Ng  are  independent,  and  requires  only 
knowledge  of  the  distribution  of  S  under  the  assumption  that  NA  and  Ng  are  independent. 
The  second  problem  is  much  more  difficult:  to  compare  S  across  two  situations  we  must 
know  the  distribution  of  S  when  the  two  point  processes  are  dependent. 
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In  a  more  general  context,  Ripley  (1976,1977)  introduced  a  measure  K,  defined 
on  an  appropriate  space,  that  summarizes  the  second-order  properties  of  the  process. 
Before  describing  this  measure,  we  need  to  state  some  assumptions  and  introduce  some 
notation.  Let  N^(s,t)  denote  the  number  of  events  of  type  i  occurring  in  the  interval 
Cs.t],  for  i  =  A,  B.  Assume  that  each  process  has  no  multiple  occurrences,  and  that 
the  intensities 

X.  =  lim  P(N. (t,t+h)  > 0}  for  i  =  A,  B  (1) 

h-K) 


are  finite.  (The  existence  of  these  limits  was  proved  by  Khintchine,  1960).  The 
X^'s  then  have  an  interpretation  as  mean  number  of  occurrences  per  unit  of  time:  for 

V*2- 

E  NjUj.tj)  =  W'l1  for  1  *  A.  8  (2) 

(this  follows  from  Dobrushin's  Lemma  and  Korolyuk's  Theorem;  see  Leadbetter,  1968). 

We  now  give  an  informal  description  of  the  measure  K,  adapted  to  the  present 
context. 

The  measure  K  is  defined  on  the  Borel  subsets  of  R,  and  for  tj  c t2,  writing 
K^tl,t2^  ^or  we  have 


K^tl’t2^  =  T~  ^^A^l’^  !a  ®  point  at  t=0) 

A 

I"  E{Ng(-t2»-t^)  |an  A  point  at  t=0}| 


C3) 


Note  that  if  NA  and  Ng  are  independent,  then 


K(tl»t2)  =  l2  ~  tl  » 


regardless  of  the  values  X,  and  X 

A  D 


Ripley  proposed  the  estimate  of  Kft^^)  given  by 


nB  nA 


K(tl,t2-)  =  N 


anb 


l  l  I(A-Be(t  t)}, 
i=l  j=l  31  1  £ 


(5)  ;  ; 

codes 

i/or 


Dlst 

1 
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where  I{*}  denotes  the  indicator  function  (actually,  the  estimate  proposed  by  Ripley 
has  an  edge  correction  for  points  near  the  boundary  of  the  period  of  observation;  this 
edge  correction  will  not  concern  us). 

Previous  work  on  the  estimator  K  is  concerned  with  spatial  processes.  The 
results  center  on  using  R  to  test  that  a  single  process  in  Poisson  (Ripley,  1977; 
Chapter  8  of  Ripley,  1981;  Silverman,  1976)  and  on  using  K  to  test  for  independence 
of  two  processes  (Lotwick  and  Silverman,  1982;  Diggle  and  Milne,  1983). 

In  this  paper  we  study  the  asymptotic  properties  of  K(tj,t2).  The  main  result 
is  that  under  certain  regularity  conditions,  as  ng  -*■ 

^  (K(tj,t2)  -  Utl,t2))  *  W(0,o2(t1,t2)),  (6) 

where 

2 

o  (t^,t2)  may  ^e  consistently  estimated  from  the  data.  (7) 

Besides  providing  the  basis  for  a  test  of  independence  between  and  Ng,  (6)  and 
(7)  enable  one  to  test  whether  or  not  K(t^,t2)  has  changed  in  the  experimental  situ¬ 
ation  described  earlier. 

The  cross-intensity  function  defined  by 

*AB(u)  =  lim  P{NA(u+t,ut-t+h1)  >  0;  Ng(t,t+h2)  >  0}  (8) 

hjjhj  +  O  1  2 

is  related  to  K  by 

K(t  t.)  =  -i- 

Vb 

Under  the  independence  hypothesis,  A^g  =  A^g.  Brillinger  (1976)  considered  the  random 
function 

nB  nA 

jJb(u)  =  l  l  HA.-Bi  e  (u-h,u+h) }  (10) 

i=l  j=l  J 

and  showed  that  under  suitable  regularity,  if  h ->■  0  and  T-*-°°  in  such  a  way  that  hT 
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T  T  T  T  T 

remains  constant,  then  for  Uy  +  u^,  |u^-u^>|  >  2h,  1  s  k  <  k' s  M,  JAfi  (u^)  are  asymptoti¬ 
cally  independent  Poisson  random  variables  with  means  2hT  »  for  ^1.  •••»  M. 

,  Jlg(u) 

Thus,  a..(u)  *  -  can  be  used  to  estimate  XAD(u)  at  a  finite  number  of  points. 

AB  2hT  AB 

*  A 

In  practice  one  would  graph  the  two  functions  X^  and  K  over  a  finite  range,  say 
C-L,L]  (i.e.  graph  K(-L,t)  for  -L  s  t  s  L) .  Although  from  a  mathematical  viewpoint 
>AB  and  K  contain  essentially  the  same  information,  the  statistical  properties  of 
their  estimates  are  quite  different:  estimation  of  X^g  is  akin  to  estimating  a  den¬ 
sity,  and  from  Brillinger's  result  the  variance  of  X^g  is  of  the  order  on  the 

other  hand,  estimation  of  K  is  akin  to  estimating  a  distribution  function,  and 

A  1  A 

from  (6),  the  variance  of  K  is  of  the  smaller  order  _.  A  graph  of  XAg  may, 

nB 

however,  indicate  features  (spikes,  location  of  maxima  and  minima,  etc.)  that  cannot 

A 

be  seen  in  the  graph  of  K.  Clearly  the  two  approaches  are  complementary. 


2.  ASYMPTOTIC  DISTRIBUTION  OF  THE  K-FUNCTION. 


UAB^tl,t2^  =  E^NA^tl't2^a  B  Point  at  t=°}' 


may  estimate  u^B^tl,t2'* 


nB  nA 


^ABftl,t2)  =  nR  "  Bi £  ^tl’t2)} 

D  1  —  1  J  =  1 


Letting 


we  note  that 


^  *  y  for  i  =  A,  B 


K^tl»t2’*  *  X.  UAB^tl,t2^'  311(1  ^tl’t2^  =  -  ^AB^tl,t2')* 

A  A\ 


To  prove  asymptotic  normality  of  K  (Theorem  2)  we  will  prove  joint  asymptotic 
normality  of  (U^gCt^ .tj) ,  X^)  .  We  will  in  fact  find  it  necessary  to  first  prove 


joint  asymptotic  normality  of  ( <^e^ta'oet^0<*  (i»e*  a 

xz 

first  term  Taylor  expansion)  applied  to  the  function  f(x,y,z)  =  —  then  yields  the 
asymptotic  normality  of  K(tj,t2).  We  also  obtain  the  joint  asymptotic  normality  of 

^AB^l*^’  *A*  *B^  by  aPP1yin8  delta-method  to  the  function  g(x,y,z)  *  (x,y/z,l/z). 

We  now  need  to  give  the  statistical  setting  of  our  asymptotic  investigation.  The 
functions  and  K(t1,t2)  involve  the  notion  of  the  Palm  measure.  That  is 

for  e  >  0,  we  consider  the  conditional  distribution  of  the  process  (NA»Ng)  given  that 
there  is  a  B  point  in  the  interval  (0,e),  and  take  the  limiting  distribution  of 
(Na,Nb)  as  e-*-0.  Intuitively,  this  corresponds  to  selecting  a  B  point  "arbitrarily", 
and  considering  the  process  with  that  point  labeled  the  origin.  This  notion  is 
discussed  for  univariate  processes  by  Leadbetter  (1972)  and  for  bivariate  point 
processes  by  Wisniewski  (1972).  We  will  assume  that  the  process  is  observed  during 
a  period  of  length  T  starting  immediately  after  the  occurrence  of  an  "arbitrary"  B 
point,  say  BQ  (thus,  we  will  be  working  with  the  Palm  measure).  This  mode  of  sampling 
is  called  semi synchronous  sampling  by  Cox  and  Lewis  (1972);  see  Wisniewski  (1972)  for 
some  fundamental  properties  related  to  it.  Also,  for  the  sake  of  convenience,  we 


will  assume  that  the  period  of  observation  ends  with  a  B  point. 


B 

Let  F 


0 


denote  the  o-field  generated  by  the  events 


{V<W  W>—*Bk  e(W  VWm); 

1  m 


+  ri ,  Bq  +  Sj)  =  hj ,  . . . ,  Na  (Bn  ♦  rn ,  Bn  +  sn)  =  hn> 


A^O  An’  “0  n'  ‘V 


for  v .  <  w .  s  0 ,  k.  =  -1,  -2,  ....  i  =  1 ,  ...,  m;  r .  <  s  .  s  0 ,  h .  =  0 ,  1,  2,  ...,  j  =  1 ,  ..., 

l  l  ’  l  *  *  ’  *  j  j  ’  j  J 

00 

n,  and  m  and  n  nonnegative  integers.  For  u>0,  let  denote  the  o-field  generated 

by  the  events 


(Bv  «  (BQ  +  Vj,  Bq  +  w^,  ....  \  e  (Bn  +  vm,  Bn  +  wJ  ; 


k  '  m’  0  nr 

m 


NA(Vrl*  B0  +  sl)=hl . NA(Vrn’  80  +  sn)=hn} 
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for  usv.<w.,  k.  =  1,  2,  ....  i»l,  ....  »;  usr.  <s.,  h.  =0,  1,  2,  ...,  j  =  1,  ....  n, 
and  m  and  n  nonnegative  integers. 

Let  a 

“n  od 

(IS) 


a(u)  =  sup{  |P(EjnE2)  -  P(EX)  P(E2)|  jE^F.^  E2  e  FB()>u} 

If  a(u)-*0  as  u-*-®,  then  the  distant  future  is  virtually  independent  of  the  past. 
We  will  actually  need  stronger  conditions  on  a(*)  . 

Let  B>0,  n  >  1,  0  <  x  <  1  be  any  constants  satisfying 


(16) 


Assumptions : 


A1  /  Ca(t)]T  tS  dt  <  ® 


A2  sup 

-«< 


A3  E 


A4  E 


up  e|[N  (j,j+l)]n  I  a  B  point  at  t' 
j  <00  l  * 


=  oj  =  M 

(  'Hi‘"  777J|  1 

|^A^tl,t2^  a  B  P°*nt  at  t=0j  <  “ 

r  4^i+ 117)*] 

L(Bi-V  J<” 

r  1  4(1+ttt) 

Lna(bo'bi)j 


A5  E 


These  assumptions  are  discussed  towards  the  end  of  this  section. 


THEOREM 


1.  Assume  A1  and  A2.  Let  ^AB^tl,t2^'  *i  ^or  *  3  E 


be 


defined  by  (11),  (12),  and  (13),  respectively. 


(i)  Under  A3,  we  have  as  nfi  -*■  « 


^  ^AB^l*^  ‘  UAB(‘tl’t2^  "*■  W(-0,Y  * 
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Furthermore,  y2(t1,t2)  can  be  consistently  estimated  from  the  data. 

(ii)  Under  A4  and  AS  we  have  as  n^  -+  ® 

^  (VXA*  VXB3"  * 

Furthermore,  A  can  be  consistently  estimated  from  the  data. 

(iii)  Under  A3  -  AS  we  have  as  n^  -*■  ® 

vSi  ^AB(-tl,t23  '  UAB*‘tl,t23,  XA-XA'  XB~XB3<  >t23  3  * 

Furthermore,  E(t^,t2)  can  be  consistently  estimated  from  the  data. 


PROOF. 


(i)  We  begin  by  showing  asymptotic  normality.  Let  U7  and  be  defined  by 


and 


Uir  *  j.L.IlVtVV’i***)>  !1B0SAjSBnB> 


Note  that 


-  l  itA  ^B.nj.B.n^}, 

J  ~ 


Vw  ■  ^  T.ur 


(17) 
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n 


B  i=l 


B 


It  is  clear  that  J  U.  -  T  Uf  =  0  (1).  Thus,  it  suffices  to  prove  the  result  with 
i-1  1  i=l  1  P 

CO 

IL's  instead  of  Uf’s  in  (17),  Observe  that  the  sequence  {ILh^  w  is  stationary, 
with  mean  (^^g(t^»t23  finite  variance  (by  A3).  The  U^'s  may  be  far  from  indepen¬ 
dent:  for  small  k,  U.  and  U.  .  may  be  nearly  identical.  If,  however,  U.  and  U-  . 
are  "nearly  independent"  for  large  k,  then  one  can  still  hope  to  have  a  Central 
Limit  Theorem  effect.  The  proof  consists  of  translating  Al,  the  mixing  condition  on 
the  point  process,  into  a  mixing  condition  on  {U^}  that  allows  the  application  of  an 
appropriate  central  limit  theorem  for  stationary  sequences. 
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Let  u(k)  be  defined  for  k  =  1,  2,  ... 


by 


u(k)  =  supdPCEjnE^  -  P(H1)P(E2)  | ;  Ej  e  o(. . .  U_j  .  UQ) ,  E2  e  o(Uk,  ,  ...)}  (18) 


(Here,  o(...U_^,Uq)  denotes  the  a-field  generated  by  { . . .  ,U_1,Uq}  ,  and  similarly  for 

o(Ujc,Uk+i> . . .)) .  The  function  y(*)  is  called  the  mixing  coefficient  of  the  sequence 

00 

(U. }.  Our  goal  is  to  prove  that  £  [y(k)]T  <  ».  It  will  be  more  convenient  however, 

®  k=l 

to  show  instead  that  £  [y(2k)]T  <  «.  The  two  conditions  are  equivalent  since  y(») 

k=l 

is  nonincreasing. 

Let  k  2  1  be  fixed,  let  E^  e  o(...,  U  k) ,  E^  e  •  •  •) »  and  consider  P(E^nE2) . 


Let 


c_k =  {B.k“ Bo Ckl/(e+1)]}  and  ck =  {Bk“Bo- tkl/(s+1)n  • 


We  may  write 

P(EinE2)  =  P{ CE1nC_k)  n  (E^)  }  ♦  PUE^)  n  (C*kuCjj))  ,  (19) 

where  c  denotes  complementation. 

Consider  the  first  term  on  the  right  side  of  (19).  For  all  large  k,  since 
B_k  -  BQ  <  -  k  implies  that  B  jc  +  t2"B0’  we  ^ave  Ei  n  C  k  e  '  Furthermore, 

E2nCk‘  FB0*tk1/(6'’1h*t1  '  11'erefore’ 

PUBjUC  ^  n  (EjOl^))  s  P(Ej)P(E2)  .a([k  tj)  .  (20) 

The  second  term  on  the  right  side  of  (19)  is  obviously  less  than  or  equal  to 
P(C*k)  ♦  P 
Consider  P(C^).  Observe  that 

P(C^)  s  P{one  of  the  intervals  (B^+j  .B^+j  +  l  ,  j  =0 , 1 , . . . ,  [k^  -  1  (21) 

has  at  least  [kB/^B+^l  points}. 


(C^) .  These  last  two  probabilities  are  dealt  with  in  the  same  way. 
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By  A2,  Chebyschev's  inequality  and  Boole’s  inequality,  the  right  side  of  (21)  is 
less  than  or  equal  to  [k1/(-S+1^]  M  [k6/('8+1-)] "n  .  Combining  this  with  (20)  and 
handling  the  opposite  inequality  in  s  similar  way,  we  obtain 

y(2k)  k  .(k^^.tj)  .  k1/(6*1)MCkB/(6*1)]"’ 


Assumption  A1  implies  that 

l  (a( k1/(l>+\)}r  <  ~  .  (22) 

k-1 

QO 

Combining  (22)  and  (16)  we  obtain  that  £  [y(2k)]T  <  »,  and  hence  that 

k=l 


l  Cy(k) ]T  <  “  .  (23) 

k=l 


Assumption  A3  implies  in  particular  that 


J1 


This,  together  with  (23)  allows  us  to  apply  Theorem  18.5.3  of  Ibragimov  and  Linnik 

2 

(1971)  to  conclude  that  for  y  (tj.t^  defined  by 

<30 

.2, 


=  Var  UQ+2  £  Cov(U(),Uh) 
h=  1 


(24) 


we  have  as  ng  -*•  ® 


\  ('5AB(tl-t2)-U*B(tX't2)> 

Consider  next  y2(t1,t2).  Let  vh  =  Co  v  (1^,1^)  and  let  vh  denote  the  sample 
covariance  at  lag  h: 

n  g  "■  h—  1 
«  1-0 


nB-h 
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where  0  =  0._(t.,t_).  Defining  the  spectral  density  of  {U.}  by 

Ad  x.  l  1 

00 

f(x)  =  7?  vo  +7.1.  \  cos  xh* 

n=  1 

2  2 

we  see  that  y  (t^,t2)  =  2uf(0) .  Thus,  to  estimate  y  (tj,t2)  we  have  available  the 
machinery  for  spectral  density  estimation  from  the  time  series  literature.  Frequently 
used  estimates  of  2irf(0)  are  of  the  form 

C0  +  2X  (25) 

h=  1 

where  s  and  c  ,  .  c  are  constants  depending  on  nD,  with  s/nD  0  as  nD  +  ®,  and 

IS  D  D  D 

C1 ~  c2  “  -  cs*  The  choice  of  s  and  (c^l^.^  is  not  discussed  here.  For  such  a 

discussion,  see  any  standard  text  on  time  series  (e.g.  Section  7.4  of  Chatfield, 

1980;  Chapter  9  of  Anderson,  1971).  Consistency  results  for  spectral  density  estimates 
have  been  established  under  conditions  on  {IK}  that  are  not  implied  by  Assumptions 
A1  -  AS  (e.g.  existence  of  all  moments  in  Brillinger,  1975;  UK}  is  a  linear  process 

as  in  Anderson,  1971  and  in  Hannan,  1970).  In  the  appendix,  it  is  shown  that  there 

2 

exists  a  consistent  sequence  of  estimators  of  y  (tpt2). 

(ii)  Let  $2  €  F,  and  let 


X. 

1 


<1  * 


«2<Bi-Bi-l> 


The  sequence  (X.}°°  is  stationary,  and  if  rr(»)  denotes  its  mixing  coefficient,  it 
1  1=-°° 

is  clear  that  (23)  holds  for  ir(*)  as  well.  This  gives  a  central  limit  theorem  for 
{X.},  and  by  the  Cramer-Wold  device  we  have  that 


VW’  nj-  E(B1 


V 


is  asymptotically  normal  with  mean  0,  and  a  covariance  matrix  which  can  be  consistently 

X  1 

estimated  as  in  Part  (i).  Applying  the  delta  method  with  the  function  f(x,y)  =  (— ,  —  ) 
we  obtain  that 
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j —  fnA  E  NA(B0,BP  nB  1  ' 

B  (T  _  E(BrB0)  *  T  '  E(BrB0) 

is  asymptotically  normal,  and  it  is  simple  to  argue  that  E  N. (Bn,B. )/E(B, -Bn)  =  X 

A  v  i  I  U  A 

and  (ECB^Bq)}-1  =  Xp  . 

(iii)  The  proof  of  (iii)  is  similar  to  that  of  (ii)  and  is  omitted. 

THEOREM  2.  Let  Kft^.t^)  and  RCt^^)  be  defined  by  (3)  and  (5),  respectively,  and 
assume  A1  -  A5.  Then,  as  n^ 

(R(t1,t2>  -  K(tx,t2)  i  N(0,o2(t1,t2)) , 

2 

where  o  (tj,t2)  can  be  consistently  estimated  from  the  data. 

Results  giving  the  asymptotic  normality  of  estimates  of  XA  and  Xg  (under  varying 
sets  of  assumptions)  already  exist  in  the  literature;  see  e.g.  Theorem  8.6  of  Daley 
and  Vere-Jones  (1972) .  It  was  necessary  to  establish  joint  asymptotic  normality  of 
XA  and  0AB(tlft2)  in  order  to  obtain  asymptotic  normality  of  K(t^,t2). 

Any  theorem  giving  asymptotic  normality  of  the  normalized  partial  sums  of  a 
stationary  sequence  {T^}  must  assume  a  moment  condition  on  T^  and  also  a  mixing 
condition  on  {T^}.  In  general,  weakening  of  the  moment  condition  must  be  compensated 
by  strengthening  of  the  mixing  condition,  and  vice  versa.  Assumptions  A3,  A4,  and 
A5  provide  moment  conditions  on  the  sequences  {IT},  {NA(B^  an<*  ^Bi.~Bi  i^» 

respectively.  Assumption  A2  insures  that  the  B  process  "moves  along"  rapidly 
enough  so  that  Al,  the  mixing  condition  imposed  on  the  point  process,  translates 
into  a  mixing  condition  for  the  sequences  {U.},  (N.(B.  .,B.))  and  { (B . - B .  .)}. 
Relationship  (16)  describes  in  a  technical  way  the  interplay  between  the  mixing  rate 
on  the  point  process  and  the  moment  condition  on  the  sequences  { ) ,  (NA(B^  ^,B^)}, 
and  {(B  -Bi 
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The  conditions  assumed  by  Brillinger  (1976)  neither  imply  nor  are  implied  by 
A1  -  AS  of  the  present  paper.  Brillinger  assumes  a  mixing  condition  on  the  bivariate 
point  process  and  also  that  the  "second  order  moments"  ^j(’)  (i,j=A,B)  exist  and 
are  continuous  (he  also  assumes  existence  and  continuity  of  the  "third  and  fourth 
order  moments";  see  equation  (2.2)  of  his  paper).  This  condition  on  XAR(«)  is  not 
satisfied  by  the  following  process:  Ng  is  a  Poisson  process,  and  is  Ng  shifted  to 
the  right  by  1  unit.  In  this  case,  This  process  does  however  satisfy 

A1-A5.  Conversely,  it  is  easy  to  find  processes  (N'A,Ng)  satisfying  all  of  Brillinger's 
conditions,  but  not  those  of  the  present  paper.  Perhaps  the  simplest  example  is  the 
following.  Let  and  Ng  be  independent,  being  a  Poisson  process,  and  being 
an  equilibrium  renewal  process  on  (-®,00)  (for  a  definition  and  a  construction  see 
pp.  517-19  of  Karlin  and  Taylor,  1975)  with  interarrival  distribution  having  a  first 
moment  but  no  second  moment.  Then  A4  is  violated,  and  it  is  not  difficult  to  check 
that  this  process  satisfies  all  of  Brillinger's  conditions. 

3.  DISCUSSION. 

The  contributions  of  this  paper  are  two-fold:  proof  of  asymptotic  normality  of 
K(ti,t2)  and  a  method  for  estimating  the  asymptotic  variance  a2(t^,t2),  enabling  the 
construction  of  asymptotic  confidence  intervals  for  K(t^,t2),  for  fixed  values  of 
tj  and  t2* 

The  function  K(*,*)  will  usually  be  of  interest  over  a  continuum  of  values,  say 
-L<tj<t2sl,  where  L  is  some  number  much  smaller  than  T.  One  can  plot  K(-L,t) 
for  -L  s  t  s  L  or,  what  is  sometimes  more  useful,  plot  K(t-y,  t+y)  for 
-L+  -  StsL-^.  Here,  d  is  some  small  number  representing  the  experimenter’s 
guess  at  the  duration  or  likely  duration  of  the  effect  of  a  B  point  on  the  A  process. 

The  function  K(t-j,  t+y)  is  identically  equal  to  d  if  NA  and  Ng  are  independent. 
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We  may  form  the  bands 


R(-L,t)  ±  z(a/2)  o(- L,t)//n^  -L  <  t  U 

and 


z(a/2) 


t  <  L 


d 


2 


where  ^Ctptj)  is  an  estimate  of  a(t1>t2),  and  z^2^  is  the  upper  j  •  100  percentile 
point  of  a  standard  normal  variable.  These  bands  of  course  are  not  simultaneous 
confidence  bands.  To  form  simultaneous  confidence  bands  one  would  need  to  carry  out 
two  distinct  steps: 


(i)  establish  weak  convergence  of  the  processes 


and 


V  (t)  =  *OK(-L,t)  -  K(-L,t) ) 
nB  B 


to  Gaussian  processes  V(t)  and  W(t),  respectively. 

(ii)  obtain  v^  and  w^° ,  the  upper  a  •  100  percentile  points  of  sup  |V(t) 

-L<tsL 

and  sup  |W(t)|,  respectively. 

-L+lst^L— 

2  2 


The  bands 


and 


K(-L,t)  ±  v(o)/C  -L  <  t  <  L 


K(t-|,  t+|)  ±w(a)/^ 


-L  4  <  t  U  •  } 
2  2 


are  then  asymptotic  simultaneous  confidence  bands. 

A  proof  of  weak  convergence  appears  extremely  difficult.  Although  desirable 
from  a  theoretical  point  of  view,  weak  convergence  is  not  useful  statistically  unless 
the  distribution  of  the  supremum  of  the  absolute  value  of  the  limiting  process  can 
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be  obtained.  In  general  this  is  a  very  difficult  problem  even  if  the  Gaussian 
process  is  stationary  (see  Cressie  and  Davis,  1981).  In  the  case  of  two  independent 
Poisson  processes  weak  convergence  of  U^C*)  an<*  Wn^(*)  can  established.  Wn^(*) 
converges  weakly  to  a  stationary  Gaussian  process,  the  distribution  of  the  supremum 
of  which  is  known.  Since  this  distribution  depends  only  on  the  two  rates  and 
it  can  be  estimated  directly  from  the  data.  These  results  will  be  reported  in  a 
future  paper. 
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APPENDIX 


In  this  appendix  it  is  shown  that  under  Assumptions  A1  -  AS,  the  asymptotic 

2 

variances  in  Theorem  1  can  be  consistently  estimated.  We  discuss  only  y  (tj^); 

the  other  asymptotic  variances  are  handled  in  the  same  way.  We  emphasize  that 

in  applications  simple  estimates  such  as  those  given  in  (25)  would  suffice. 

Consider  the  sequence  {U^}.  For  an  arbitrary  integer  n,  let 

n 

l  U. 

Y.  =  .  (26) 

/n 

It  is  well-known  (and  easy  to  see)  that 

Var  Y.  =  I  (l-~)  v„  =  y2(t.,t_)  +  e(n),  where  e(n)  +  0  as  n  +  ».  (27) 

1  j£|<n  n  l  12 

If  Yj,  Y2,  ....  Yh  are  i.i.d.,  then 

“h  -  “2 

l  (Y.-EY.T 

E  i zl - Y2(tj,t2)  =  e2(n)  +  i  Var(Y1-EY1)2  .  (28) 

h 

By  A3,  VarfY^-EY^2  is  finite. 

nB 

The  idea  is  to  divide  the  sequence  into  h  blocks,  each  of  size  n, 

with  the  blocks  separated  by  a  distance  f,  and  to  let  Y^  be  the  normalized  sum 
of  the  U's  in  the  jth  block,  as  in  (26),  for  j  »  1,  2,  ...,  h;  see  the  diagram. 


Diagram  to  describe 
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Then 


-T  "  IT  (V*)2*  where  ?  =  K"  (29) 

can  be  used  to  estimate  y  The  proof  consists  basically  of  making  the 

following  heuristics  rigorous.  If  f  is  large,  then  by  the  mixing  condition  on 

the  sequence  (U^)  (see  (23))  the  Y^'s  will  be  nearly  independent,  so  that  equation 

(28)  will  be  approximately  true;  if  n  is  large,  then  e(n)  in  (27)  and  (28)  will 

be  small,  and  if  h  is  large  then  the  second  term  on  the  right  side  of  (28)  will 

2 

be  small  (note  however  that  Var(Y^-EY^)  depends  on  n) .  Also,  the  effect  of 

substituting  Y  for  EY^  in  the  left  side  of  (28)  should  be  negligible.  Thus,  T 
2 

converges  to  y  (t^,t2)  in  probability. 

The  numbers  n  =  n(ng) ,  h  =  h(ng)  and  f  =  f(n^)  of  course  will  depend  on  n^, 
but  in  a  way  that  is  to  be  determined  later.  The  dependence  on  n^  will  be 
suppressed  in  the  notation  for  convenience. 

We  begin  by  decomposing  T  given  by  (29)  as 

t4  1  (Y.-EY.)2  -  (Y-EY  ) 2  .  (30) 

h  i=l  1  1  1 

2 

We  will  show  that  the  first  term  in  the  right  side  of  (30)  converges  to  y  (t^t^ 

in  quadratic  mean  (hence  in  probability),  and  that  the  second  term  in  the  right 

side  of  (30)  converges  to  0  in  mean  (hence  in  probability) . 

Considering  first  \  (Y.-EY. )2,  we  have  by  (27) 
h  i=l  1  1 

-  h  2  h 
E  -  I  (Y.-EY.)2- y2(t1,t2)  =  e2(n)  +  Var j  £  \  (Yi-EY±) 2  l  .  (31) 

_h  i*l  J  i=l 

Writing 

V  =  ivar(Y.-EY  )2  and  C  =  4ll  Cov((Y.-EY  )2,  (Y-EY.)2)  (32) 

h  1  1  h  i*j  11  J  J 
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we  have 


Var 


ir  I  <VEV: 

1  =  1 


=  V  +  c  , 


(33) 


and  we  examine  V  and  C  separately,  beginning  with  V. 

Bounding  the  variance  of  a  random  variable  by  its  second  moment  and  using 
(26) ,  we  may  write 


v  s  — -y  e 

hn2 


it 

l  U. 
i 


Li=1 


(34) 


The  Minkowski  inequality  gives 


V  s  -L.  n4Eu5  =  -  0(n2) 
hn2  1  h 


(35) 


by  A3. 


We  now  consider  C  in  (33) .  If  i  <  j  are  fixed  and  Us  and  Ut  are  summands  of 
/n  Y^  and  i/n  Y^,  respectively,  then  t  2  s  +  f .  Therefore,  by  Theorem  17.2.2  of 
Ibragimov  and  linnik  (1971) 


|Cov((Yi-EY.)2,  (Y.-EY.)2)|  S  Cu(f)]T  |4  +  6 


1-T 

/  I  l4(1  +  T^T)\  2 

{e|VeyiI  ) 


,  (36) 


with  u(*)  defined  by  (18).  Now  by  the  Minkowski  inequality. 


E|Y1-EY1 


..  .. 


s  n 


ElVERi' 


(37) 


and  therefore,  (36)  can  be  rewritten  as 

|Cov((Yi-EYi)2,  (Yj-EYj) 2)  |  =  Cu(f)  ]T0(n) . 
Since  the  sequence  { Cvi(f)  3T) £_ ^  is  nonincreasing,  we  may  write 


(38) 


17 


(39) 


Such  sequences  are  very  easy  to  construct.  Then  the  right  sides  of  equations 
(35),  (40)  and  (42)  converge  to  0  as  nfl  +  ®.  This  completes  the  proof. 
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